Spin correlations and exchange in square lattice frustrated ferromagnets 
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The J1-J2 model on a square lattice exhibits a rich variety of different forms of magnetic order 
that depend sensitively on the ratio of exchange constants J2/ Ji . We use bulk magnetometry and 
polarized neutron scattering to determine Ji and J2 unambiguously for two materials in a new 
family of vanadium phosphates, Pb2VO(P04)2 and SrZnVO(P04)2, and we find that they have 
ferromagnetic Ji. The ordered moment in the collinear antiferromagnetic ground state is reduced, 
and the diffuse magnetic scattering is enhanced, as the predicted bond-nematic region of the phase 
diagram is approached. 

PACS numbers: 75.50.y; 75.40.Cx; 75.40.Gb; 75.30.Et; 75.25.+Z 
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The square-lattice S = 1/2 Heisenberg antiferro- 
magnet with nearest-neighbour (NN) exchange constant 
Ji and next-nearest-neighbour (NNN) exchange J2, has 
long served as a paradigm for the two-dimensional frus- 
trated magnetism (for a review see Ref. [1]). It is de- 
scribed by the following Hamiltonian: 



— Ji y ^ Si • S j -\- J2 ^ ^ Si ■ 



(1) 



<ij>i 



<i,k>2 



where < i, j >i and < i^k >2 refer to pairs of NN and 
NNN respectively. The so-called Ji — J2 model first 
came to prominence due to its relevance to the high- 
temperature superconducting cuprates, and the recent 
discovery of high- Tc superconductivity in pnictides [l, Q 
has increased interest in this model, since it may play a 
key role for the magnetism on the square Fe sublattice [J] . 

The frustrated square lattice model is characterized by 
the frustration ratio a = J2/J1, and the energy scale is 
given by Jc = + J| . When the exchange constants 
are antiferromagnetic (AF) and a ^ 0.5 the ground state 
is believed to be a valence bond solid (VBS) in which 
spins form tightly bound singlets on nearest-nieghbour 
bonds On the other hand, recent calculations for 
ferromagnetic (FM) Ji and AF J2 predict d-wave bond- 
nematic order for a ^ —0.5 0, 0- The calculated phase 
diagram is shown schematically in Fig. [TJ 

There are now several experimental realizations of 
spin- 1/2 magnets on a square-lattice where both Jiand 
J2 interactions play an important role. The first ex- 
amples, Li2VOX04 (X = Si, Ge), are collinear anti- 
ferromagnet s (CAF) with AF exchange J2 ^ Ji [1, 
i, [13, O, EH, H Q whereas VOM0O4 and 
PbVOs are Neel antiferromagnets (NAF) with AF ex- 
change J2 <C Ji l3, [l^- In contrast, the new family of 
vanadium phosphates AA'VO(P04)2 (A,A' = Pb, Zn, Sr, 



Ba. Cd) [Ti,0,[2ll,[22[ and (CuX)LaNb207 (X = CI, Br) 
[23I . [2^ have been argued to be frustrated ferromagnets 
with FM Ji and AF J2. 

In order to explore this phase diagram, accurate con- 
trol of a is required. Until now it has proved difficult to 
determine the exchange constants accurately using bulk 
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FIG. 1: Calculated phase diagram for the general Ji — 
J2 model, locating Pb2VO(P04)2 and SrZnVO(P04)2. The 
values of a in parentheses are thephase boundaries deter- 
mined by exact diagonalization 0, At these values, zero- 
point fluctuations destroy the relevant order parameter. The 
yellow area corresponds to the valence bond solid (VBS) 
regime, the purple zone to the bond-nematic (d-BN) region 
0? 0, [2^. The straight lines are obtained from the Curie- 
Weiss temperatures, and the circles are obtained from the 
maxima in susceptibility, as described in the text. The filled 
circles are the solutions selected using polarized neutrons. 
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thermodynamic properties, and this is reflected by the 
wide variation reported for a for individual compounds 
in the literature. The determination of a from thermo- 
dynamic quantities alone is undermined by the fact that 
J I + J 2 and I Ji + J2 1 , but not Ji — J 2 are well determined. 

Here we present combined bulk magnetization and po- 
larized neutron scattering studies of Pb2VO(P04)2 and 
SrZnVO(P04)2. The determination of the magnetic 
structures in the ground state, and the comparison of the 
diffuse neutron scattering intensity in the paramagnetic 
phase with high-temperature series expansions (HTSE) 
provides an unambiguous determination of the exchange 
constants Ji and J2 . We are also able to explore the con- 
sequences for the spin correlations of the approach to a 
quantum disordered region of the phase diagram. 

Large powder samples were synthesized at the Max- 
Planck Institute for Chemical Physics of Solids in Dres- 
den by the solid state reaction method, in order to 
avoid multi-phase composition for these compounds. The 
quasielastic neutron experiments were performed on the 
D7 diffuse scattering spectrometer at the Institut Laue- 
Langevin (ILL) in Grenoble, with a fixed incident wave- 
length of A = 3.1 A corresponding to an energy window of 
8.5 meV. XYZ polarization analysis was employed to sep- 
arate magnetic signal from the coherent structural and 
incoherent scattering [26|]. Polarized inelastic neutron 
scattering measurements were performed using the IN20 
triple- axis spectrometer of the ILL, in order to confirm 
the energy scale of the system. Scans of energy transfer 
were performed at various reciprocal lattice points, Q, 
and temperatures, T, and in all cases the excitations were 
found to be within ±4 meV. Thus the measurements on 
D7 integrate over the full spectral energy range to yield 
the structure factor, S{Q). The bulk magnetisations of 
Pb2VO(P04)2 and SrZnVO(P04)2 were measured using 
a SQUID magnetometer at the University of Liverpool. 

Pb2VO(P04)2 crystallizes in the monoclinic system 
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FIG. 2: The structure of Pb2VO(P04)2 contains corrugated 
layers formed by VO5 square pyramids (blue), oriented in al- 
ternating directions in a chessboard fashion. The pyramids 
are connected by tetrahedral PO4 groups (yellow). The lay- 
ers are separated by the Pb atoms (green) and isolated PO4 
tetrahedra. Projections shown are for a) ab and b) be planes. 



with spacegroup P2i/a and lattice parameters a = 
8.747(4) A, 6 = 9.016(5) A, c = 9.863(9) A and 
P = 100.96(4)° [iO, [23. Two projections of the 
Pb2VO(P04)2 structure are shown in Fig. [2l The struc- 
ture contains corrugated layers formed by VO5 square 
pyramids oriented in alternating directions in a chess- 
board fashion. The pyramids are connected by tetra- 
hedral PO4 groups. The layers are buckled along the 
b-axis. The square bases of the similarly oriented pyra- 
mids are located approximately at the same level. The 
magnetic [VOPO4] layers are well separated (~10A) 
by non-magnetic Pb atoms and isolated PO4 tetrahe- 
dra, indicating a quasi-2D character for the system. 
SrZnVO(P04)2 shares many characteristics with its rel- 
ative Pb2VO(P04)2. Their structural relationship is 
the replacement of the two Pb^+ cations by Sr^+ and 
Zn^+. Again, it forms a two-dimensional square lat- 
tice of ions with competing interactions leading to 
frustration. The profound structural similarities with 
Pb2VO(P04)2 attracted our interest since small changes 
introduced by the cationic substitution can tune the 
magnetic exchanges Ji and J2 and thus the frustra- 
tion ratio. SrZnVO(P04)2 crystallises in the orthorhom- 
bic system with spacegroup Pbca and lattice parameters 
a = 9.0660 A, 6 = 9.0117 A and c = 17.5130 A 
The main difference between the two compounds' space 
groups arises from the way the VOPO4 layers are stacked 
along the c-axis. 

The polycrystalline samples were measured under an 
applied field of 1000 Oe, giving the results of Fig. [3l 
For T>50 K, x(T) is characterised by a paramagnetic 
Curie- Weiss (GW) behaviour. Below that, x(T) passes 
through a broad maximum that is characteristic of a low- 
dimensional magnet, peaked at TMax ~ 8.5 K and 6.3 K 
for Pb2VO(P04)2 and SrZnVO(P04)2 respectively. The 
good quality of the samples is evident by the absence 
of any trace of impurity tails at the lowest attainable 
temperatures, which is typical of paramagnetic foreign 
phases. Kinks in x(T) indicate the onset of magnetic or- 
dering at the Neel temperatures Tn ~ 3.7 K and 2.6 
K for Pb2VO(P04)2 and SrZnVO(P04)2 respectively. 
The fits of the Gurie- Weiss law to the high-temperature 
data (T>50 K) are shown as insets in Fig. [3l The ef- 
fective moments /ie// = (1.72 ± 0.02) are consistent 
with the predicted value for S = 1/2, /ig// = 1-73. 
The extrapolated intercepts with the temperature axis 
give Ocw = Ji^ J2 = (4.5 ± 1) K and (1.5±1) K for 
Pb2VO(P04)2 and SrZnVO(P04)2 respectively. The val- 
ues of ^cw are rather small compared to the positions of 
the maxima in the susceptibilities TMax, which allow the 
overall energy scales of the systems Jc to be estimated. 
This suggests that Ji and J2 have opposite signs, provid- 
ing evidence of mixed FM and AF exchange couplings . 
The conclusions from the susceptibility measurements for 
the two compounds are summarized in the phase diagram 
of Fig. [H The observed value of Oqw defines a straight 
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FIG. 3: Magnetic response of poly crystalline (a) 
Pb2VO(P04)2 and (b) SrZnVO(P04)2 as measured at 
1000 Oe. The broad peaks of x(T) are characteristic for low 
dimensional spin systems. The magnetic phase transitions at 
Tn manifest as a change of slope in x(T), while TMax indicate 
the energy scale for these systems. Insets show inverse sus- 
ceptibilities for T > 50 K, together with the Curie- Weiss law 
fits (solid lines), giving ^cw = (4.5±1) K and ^cw = (1.5itl) 
K for Pb2VO(P04)2 and SrZnVO(P04)2 respectively. 



line on this plot, while the relation between TMax and Jc, 
which has been calculated in Ref. 0, enables us to get a 
pair of solutions on this line. Thus for both compounds, 
there is one possible solution in the CAF and another in 
the NAF region of the phase diagram. 

The neutron scattering experiment complemented 
these bulk property results, yielding a unique solution 
in each case. All of the scattering in Fig. [H is mag- 
netic, and it has been placed on an absolute scale by 
reference to a standard vanadium sample. At tempera- 
tures below the Neel point, the magnetic scattering con- 
sists of sharp magnetic Bragg reflections on top of a dif- 
fuse magnetic background. Data for powdered samples 
of Pb2VO(P04)2, SrZnVO(P04)2, as well as a control 
sample of Li2VOSi04 were collected at 1.5 K and were all 
found to have a CAF structure. This immediately selects 



only one of the two solutions consistent with the bulk 
magnetization data. Fig. [4^ shows the magnetic ground 
state diffraction pattern of Pb2VO(P04)2, modelled by 
both a CAF structure with an ordering wavevector of 
Q = (tt, 0) (solid line) and a Neel structure for compari- 
son (dashed line that fails to reproduce the data). The or- 
dered moments in the CAF structures are (0.42±0.04/iB) 
for SrZnVO(P04)2, (0.50^0.04//^) for Pb2VO(P04)2, 
and (0.55±0.04/iB) for Li2VOSi04. As the ordered mo- 
ment in the ground state decreases for this series of com- 
pounds, the integrated intensity under the diffuse com- 
ponent was found to increase. 

Exchange constants are usually determined via the ex- 
citations from the ordered phase. However, the large 
single-crystal samples required for inelastic neutron scat- 
tering were not available for any of these compounds. 
The novelty in our approach was to model S{Q) in poly- 
crystalline materials by using HTSE of the static sus- 
ceptibility. In the determination of the NN and NNN 
exchange constants we used the 1*^ order HTSE given 
by: 



S{Q)^SiS+l) 



S{S + l) J{Q) 
3 ' kT 



■F\Q) (2) 



where S is the spin of the magnetic ions, F'^{Q) the mag- 
netic form factor of the individual magnetic ions and 
J{Q) is the Fourier transform of the exchange interac- 
tions Fig. [4)3 shows the diffuse neutron scattering 
intensity at T = 20 K, a temperature well into the para- 
magnetic phase for this compound. Although magnetic 
long range order has been destroyed at this temperature, 
the kT term has not completely washed out all magnetic 
interactions. The presence of oscillations in the data 
indicates the existence of short-range spin correlations. 
Note that we have measured and fixed F'^{Q) directly 
at room temperature, where no oscillations are present 
and these systems are ideal paramagnets. The calculated 
scattering from the two pairs of J's in Fig. [T] are shown 
in Fig. [4)3. Our results again select the unique solution 
in the CAF region, with FM Ji = (-3.2 ± 1) K and 
AF J2 = (7.7 ± 1) K. Similarly, for SrZnVO(P04)2, the 
combined neutron and suscebtibility measurements yield 

= (-4.6 ±1) K (FM) and J2 = (6.1 ±1) K (AF). Note 
that a similar analysis for Li2VOSi04 was not possible 
due to absorption of neutrons by the sample. 

The ordered moments in the ground states are smaller 
than the nominal value for 5* = 1/2 systems due to 
quantum fluctuations. The value for Li2VOSi04 is 
consistent with Ref. [lH. However, the values 
for Pb2VO(P04)2 and SrZnVO(P04)2 are significantly 
lower, indicating an increase in quantum disorder. This 
is entirely consistent with the fact that values of Ji and 
J2 determined in this study place them successively closer 
to the disordered region of the phase diagram. 

In summary, we have unambiguously determined the 
exchange constants that are key to understanding a new 
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ther experimental study, for example in high magnetic 
fields Sfl, mi , and the polarized neutron scattering tech- 
niques described here open new vistas for the exploration 
of these model magnetic systems. 
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FIG. 4: (a) Magnetic ground state diffraction pattern of 
Pb2VO(P04)2 measured at 1.5 K. The ordered component 
of the ions' moments gives rise to sharp magnetic Bragg re- 
flections, while the disordered component gives a diffuse back- 
ground. Fits are a CAP model with Q = (vr, 0) (solid line) as 
well as a Neel model (dashed hne, for comparison only). Part 
(b) shows short-range spin correlations in the paramagnetic 
phase of Pb2VO(P04)2. The solid curve is a first-order HTSE 
calculation for the CAP solution with PM Ji = (-3.2 ± 1) K 
and AP J2 = (7.7=bl) K, and the dashed is the NAP solution. 
The data clearly favor the CAP solution. 



region of the Ji — J2 phase diagram with FM NN in- 
teractions and AF NNN interactions. We are able to 
tie down the exchange constants using novel polarized 
neutron scattering techniques from powdered samples. 
The reduced ordered moment and increased diffuse mag- 
netic signal indicate the approach to a quantum dis- 
ordered spin-nematic region. It should be possible to 
tune the exchange to be closer to the quantum disor- 
dered region by other chemical substitutions, or by the 
application of pressure or magnetic field. These low en- 
ergy scale systems are particularly well suited for fur- 
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